Slave rotor theory of antiferromagnetic Hubbard model 
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The slave-rotor mean-field theory of Florens and Georges is generalized to the antiferromagnetic 
phase of the Hubbard model. An effective action consisting of a spin rotor and a fermion is derived 
and the corresponding saddle-point action is analyzed. Zero-temperature phase diagram of the 
antiferromagnetic Hubbard model is presented. While the magnetic phase persists for all values 
of the Hubbard interaction U, the single-particle spectral function exhibits a crossover into an 
incoherent phase when the magnetic moment m (and the corresponding U values) lies within a 
certain window m c < m < 1 — m c , indicating a possible deviation from the Hartree-Fock theory. 
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PACS numbers: 

The Hubbard model has received a lot of attention 
theoretically as a prototypical model for strong elec- 
tron correlation in low dimensions 1 . In particular, 
much has been learned about the quantum phases and 
the zero-temperature transitions between them in this 
model. Brinkman-Rice (BR) theory 2 predicts a metal- 
to-paramagnetic insulator transition at a finite Hubbard 
interaction strength U, and further refinement by dynam- 
ical mean-field theory (DMFT) in recent years supports 
the original BR picture 3 . 

Recently Florens and Georges (FG) 4 introduced an in- 
teresting re-formulation of the Hubbard model. In their 
slave-rotor (SR) theory, in similar spirit to the slave- 
boson representation, an electron operator is decomposed 
as the product of a fermion and a U(l) "slave-rotor" 
operator. Analysis of the effective mean-field action re- 
vealed that a quantum phase transition takes place be- 
tween metallic and paramagnetic insulating states as U is 
increased beyond a threshold value U c . A good qualita- 
tive agreement between FG's slave-rotor mean-field the- 
ory (SRMFT) with the DMFT predictions was achieved 
while avoiding the use of heavy numerical machinery of 
the latter method. More recently, SRMFT was employed 
in the understanding of frustrated Hubbard model on the 
triangular lattice 5 . 

Largely ignored in the above-mentioned theories 2-6 is 
the spin degrees of freedom responsible for the antifcrro- 
magnetism. The nesting of the half-filled Fermi surface 
and the onset of spin density wave are usually treated 
in the Hartree-Fock (HF) theory while the strong effects 
of Gutzwiller projection on the HF ground state are ig- 
nored. A notable exception in the efforts to go beyond the 
Hartree-Fock picture to understand the magnetic phase 
is given by the four-boson formulation of the Hubbard 
model by Kotliar and Ruckenstein (KR) 7 . In KR's theory 
strong on-site correlation effects as well as the magnetic 
order were treated at the mean-field level. Later quan- 
tum Monte Carlo study confirmed much of the mean-field 



conclusions of KR 8 . Both references 7,8 focused on the 
overall phase diagram, and the nature of the quasiparti- 
cle states in the magnetic phase was not thoroughly dis- 
cussed. A more recent study on this subject 9 concluded 
that at strictly zero temperature, as in the Hartree-Fock 
(HF) picture, the quasiparticles in the antiferromagnetic 
phase of the Hubbard model remain coherent. 

Given the new machinery of SRMFT, we feel that it 
is worthwhile to re-visit this issue in more detail. In this 
paper, we present a natural extension of the original 
SRMFT theory that allows one to treat the magnetic as 
well as the non-magnetic phase of the Hubbard model. 
The saddle-point analysis of the effective action for the 
half-filled model gives two phases, characterized by the 
coherence/incoherence of the quasiparticles. Further 
consideration of gauge fluctuation renders the phase 
transition into a crossover. The magnetic ordering 
persists for all values of U as in the HF theory. 

We start with the Hubbard model 

H = -t^ c] a C ta - (J, 4a°i<r + U J2 4r c *T c li c a C 1 ) 
ija ia i 

defined on the two-dimensional square lattice. The 
Hubbard-C/ term can be decomposed into charge and spin 
channels in standard fashion 
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with riiu = c irT Ci a . The latter term in Eq. (2) was ignored 
in the study of spinless states by FG. 

Florens and Georges introduced a decomposition of the 
electron operator as Ci a = e~ l6i fi a where the rotor vari- 



able 
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serves to keep track of the charge number at 



each site, as the creation or annihilation of an electron is 
accompanied by the phase change e ±lSi . Here we propose 
that a second rotor variable e ^* can be introduced for 
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the bookkeeping of spin numbers. The representation of 
the electron operator we propose is 



„. — p -i6i-i<J<t>i c 
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The factor a = ±1 is used to distinguish the creation of 
up and down spins along the z-axis. Although we are 
focusing on the U(l) case here, a fully SU(2)-invariant 
representation of the spin sector is also possible. 

After the substitution made in Eq. (3), the local charge 

4a c i° and tne local z-spin vela ™ appearing in 
the decomposition of the Hubbard term are replaced 
by the conjugate operators of (6i,4>i), which we denote 
(Li, Mi) = (-id/d6i,-id/d(f>i). The Hubbard Hamilto- 
nian takes on the expression 
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supplemented by the constraints Li = fj a fia and 

Mi — ^2 a °~fiafi<7 at every site 4 . The corresponding ac- 
tion can also be derived straightforwardly: 

ia ija 
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(5) 

A pair of Lagrange multiplcr fields (A^, 77^) has been in- 
troduced to impose the constraints. The original theory 
of FG is based on this Lagrangian, without the terms 
pertaining to the spin decomposition such as M i; </>j, and 
r/i. Our formulation thus generalizes the scheme of FG 
in a natural way to include magnetic order. 

To proceed further, we ignore the terms pertaining to 
the charge sector as they are already discussed by FG. 
Leaving out the terms containing (9i, Li, Xi), the next 
step is to integrate out Mi from the action. Because of 
the negative norm in M? (third line of Eq. (5)), we use 
the following identity of the Gaussian integration 
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to first re-write the action (5) with the positive norm for 
Mf. The integration over Mi of the modified action can 
then proceed to give 
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The final manipulation involves the shift 77, — > r]i + id T 
that gives 
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This concludes the formal derivation of the effective ac- 
tion for the spin-ful Hubbard model within the slave- 
rotor framework. Without the phase fluctuations in Eq. 
(8) this effective action at the saddle-point level is equiv- 
alent to the HF theory of the Hubbard model. 

To explore the consequences of phase fluctuations in 
the action (8), we replace rji by its saddle-point value rj^ 
given by 



u 



(9) 



which is identical to the saddle-point equation in the ab- 
sence of the phase fluctuation. We will write r/i to refer 
to this average value in the rest of the paper. 

One can decompose the hopping term in Eq. (8) by 
a pair of Hubbard-Stratonovich (HS) fields (ay,/?^) 5 , 
resulting in the effective Lagrangian 



L = L + Lf + Lfa 

L = tJ2( a ijPij + <*ij0ij) + JJ E ^' 
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At the saddle-point level, the HS parameters take on the 
average values = (/* T / jT + fjjn), Pij = (e"*e*). 
Assuming real and uniform mean-field solutions a^- = a, 
Pij = P, and staggered effective magnetic fields r\i = 
(— I) 1 ?/, = (-l)'m, we arrive at the mean-field effec- 
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tive action: 



Lmf — Lf + L<p, 

ia ija 
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We assumed the case of no broken time-reversal sym- 
metry, (d T fa) = 0. The fermion action is in the stan- 
dard HF form except for the renormalization of the band- 
width t — ► f3t. At zero temperature the fermion sector 
thus always remains in the magnetic phase with a gap 
to quasiparticlc excitations. The boson action, on the 
other hand, is the standard XY action modified by the 
Berry phase term iJ2i m id T fa- We analyze each of the 
mean- field actions derived in Eq. (11), beginning with 
the fermion sector. 

In analyzing the fermionic mean-field action we confine 
our attention to half-filling for which /i = 0. The mean- 
field conditions for 77 and a at T = read 
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with Q — (tt,7t). In writing down the Fourier form we 
assumed a uniform mean-field solution i(qi) = q. The 
bosonic fc-sum is also over the reduced Brillouin zone. 
The boson part can be diagonalized using a pair of oper- 
ators (7ifc„,72fc„) related to {Y kv ,Y k+Qv ) by 

Y kv = -^=(cosh(9 fe - sinh0fc)(7i fel , + 7 2 fc„) 
Y k+ Q v = -^(cosh0 fe + sinh0 fc )(7i fel ,-72fc„). (14) 

After taking cosh2# fe = q/E b k , sinh2# fe = ae k /E b k , and 
E k = \/l 2 ~ ( ae k) 2 , one gets 




Here e k = —tJ2- ei e lk ^ rj ^ r ^ (j runs over all nearest 
neighbors of 1) is the bare band in the absence of exchange 
splitting introduced by non-zero 77, and z is the lattice co- 
ordination number. The fc-sum in both equations is over 
the reduced Brillouin zone and E k — y/ ((3e k ) 2 + rf is the 
fermionic energy with a gap set by -q. 

The bosonic action derived in Eq. (11) is invariant 
under nn-raorirnl-m followed by a shift of one 
lattice spacing. The Berry phase term vanishes for m = 
and m — I leaving only the classical XY model with an 
ordered phase fa — <po at zero temperature. One might 
thus expect that small deviations such asrawO and m w 
1 still gives the ordered phase while m w 0.5 introduces 
enough perturbation to induce phase disordering. We 
present a calculation which confirms this expectation. 

Following the treatment of FG 4 we introduce the uni- 
modular field Y> to make the replacement e 1 ^* = Y{. An 
additional Lagrange multiplier qi is required to impose 
the uni-modular constraint. This extension allows us to 
examine the bosonic action at the saddle-point level. The 
effective boson Lagrangian written in terms of Yi reads 



L Y = -im^v^l^ikv - 72^72^) 
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kv 

The boson spectrum E\ is gapped if q — aD > while 
q — aD = leads to the condensation of Y. Here D = zt 
is half the bare bandwidth. Two additional relations are 
obtained at T — from the constraints (Y*Yi) = 1 and 



Solving Eqs. (12) and (16) simultaneously renders the 
self-consistent parameters (a,(3,m,q) for a given U/D. 
The Bosc condensation occurs when q — aD, which gives 
m c w 0.69. 

To get a better idea on the analytical structure of the 
set of self-consistent equations obtained above, we first 
re-write Eqs. (12) and (16) as the integration over the 
energy with a certain density of states D(e), and approx- 
imate it with a constant value, D(e) = l/(2D). The 
mean-field equations are then given by 



4 




FIG. 1: (color online) Plot of mean- field parameters 
(a, (3,m,q) for U/D obtained from solving Eq. (17). q equals 
aD at (U/D) c ~ 0.81 indicated by a dashed vertical bar. 
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The Bose condensation occurs at m c = 7r/4 ps 0.8 
((U/D) c ps 0.81), similar to the earlier estimate. The 
plot of all four self-consistent parameters with varying 
U/D using Eq. (17) is shown in Fig. 1. 

The set of mean- field equations (17) yields to = 1/2, 
not to = as in the HF theory, when U/D — > 0. We 
believe this is an artifact of the replacement e l ^ i — > Y{. 
The boson action in Eq. (11), being equivalent to the 
XY model, must show ordered phase at m — and to = 
1, and the only chance of a disordered phase occurs for 
intermediate to values. Combining this general argument 
with the mean- field analysis in terms of Y, we conclude 
that the incoherent phase exists within the window 1 — 
to c < m < m c with m c ss 0.8. In this intermediate phase 
the boson gap At, = q 2 — (aD) 2 remains finite. Due to 
the gap, the spectral function of the composite electron 
operator is incoherent 10 . The region outside this value 



gives Bose condensation, or the coherent quasiparticles. 
Meanwhile the magnetic sector remains ordered for all 
U/D as in the HF theory. For comparison we recall that 
in the paramagnetic Hubbard model 4 , SRMFT gave Bose 
condensation for small- U: < U/D < (U/D) c . 

So far we performed saddle-point analysis and obtained 
a mean-field picture showing a second order phase transi- 
tion (with the boson gap as the order parameter) for the 
spin phase field fa in the intermediate values oiU/D. It 
is natural to ask the stability of the mean-field picture 
against the gauge field ay that appears in the phase fluc- 
tuations of the hopping order parameters, ay = ae m,j 
and Pij = f3e laij , where a and (3 are the mean field val- 
ues obtained before. It should be noted that the U(l) 
spin-gauge field ay is compact, thus allowing instanton 
excitations 11 . From the seminal work of Fradkin and 
Shenker 12 we know that there can be no phase transi- 
tion between the Higgs and confinement phases due to 
instanton proliferation, and only a crossover behavior is 
expected. In the present problem the phase-coherent 
state corresponds to the Higgs phase while the phase- 
incoherent state coincides with the confinement phase. 
Applying Fradkin and Shenker's result to the present 
problem, we conclude that the second order phase tran- 
sition turns into a crossover between the coherent and 
incoherent phases. The magnetic order parameter, be- 
ing a gauge-invariant quantity, remains unaffected by the 
gauge fluctuation. 

In summary we have developed an extension of the 
slave-rotor theory of Florens and Georges to the magnet- 
ically ordered phase by introducing a second rotor vari- 
able pertaining to the spin degrees of freedom. On per- 
forming a saddle point analysis we uncover an incoherent- 
to-coherent crossover within the half-filled antiferromag- 
netic phase of the Hubbard model at zero temperature. 
The incoherent phase exists in the intermediate values of 
U/D between weak and strong coupling limits. Given the 
common conception that the magnetically ordered phase 
of the Hubbard model at T = is well understood within 
the HF theory, the possibility we suggest in this paper is 
tantalizing. The formalism developed in this work may 
also be of use for understanding other exotic magnetic 
systems. 
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